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Abstract 

In this paper, the general formulation for inextensible flows of curves 
on oriented surface in R 3 is investigated. The necessary and sufficient con- 
ditions for inextensible curve flow lying an oriented surface are expressed 
as a partial differential equation involving the geodesic curvature and the 
geodesic torsion. Moreover, some special cases of inextensible curves on 
oriented surface are given. 
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1 Introduction 

It is well known that many nonlinear phenomena in physics, chemistry and 
biology are described by dynamics of shapes, such as curves and surfaces. The 
evolution of curve and surface has significant applications in computer vision 
and image processing. The time evolution of a curve or surface generated by 
its corresponding flow in IR 3 -for this reason we shall also refer to curve and 
surface evolutions as flows throughout this article- is said to be inextensible 
if, in the former case, its arclength is preserved, and in the latter case, if its 
intrinsic curvature is preserved. Physically, the inextensible curve flows give 
rise to motions in which no strain energy is induced. The swinging motion of 
a cord of fixed length, for example, or of a piece of paper carried by the wind, 
can be described by inextensible curve and surface flows. Such motions arise 
quite naturally in a wide range of the physical applications. For example, both 
Chirikjian and Burdick [1] and Mochiyama et al. [2] study the shape control of 
hyper-redundant, or snake-like robots. 



The inextensible curve and surface flows also arise in the context of many prob- 
lems in computer vision [3], [4], computer animation [5] and even structural 
mechanics [5] . There have been a lot of studies in the literature on plane curve 
flows, particularly on evolving curves in the direction of their curvature vector 
field (referred to by various names such as " curve shortening" , flow by curva- 
ture" and "heat flow"). Particularly relevant to this paper are the methods 
developed by Gage and Hamilton 7 and Grayson [8] for studying the shrinking 
of closed plane curves to circle via heat equation. 

The distinction between heat flows and inextensible flows of planar curves were 
elaborated in detail, and some examples of the latter were given by 9.. Also, a 
general formulation for inextensible flows of curves and developable surfaces in 
R 3 are exposed by [ID) . 

In this paper, we develop the general formulation for inextensible flows of curves 
according to Darboux frame in R 3 . Necessary and sufficient conditions for an 
inextensible curve flow are expressed as a partial differential equation involving 
the geodesic curvature and geodesic torsion. 



2 Preliminaries 



Let S be an oriented surface in three-dimensional Euclidean space E 3 and a (s) 
be a curve lying on the surface S. Suppose that the curve a (s) is spatial then 

there exists the Frenet frame { f, iv\ ^} at each points of the curve where f 

is unit tangent vector, is principal normal vector and B is binormal vector, 
respectively. The Frenet equation of the curve a (s) is given by 

T 1 = kJ? 

N' = -nf + tB~ 

B' = 

where k and r are curvature and torsion of the curve a (s), respectively. 
Since the curve a (s) lies on the surface S there exists another frame of the curve 

a (s) which is called Darboux frame and denoted by ~$ i ~t%\- In this frame 

is the unit tangent of the curve, ~rt is the unit normal of the surface S and 
is a unit vector given by ~~c{ = it x T . Since the unit tangent T is common 
element of both Frenet frame and Darboux frame, the vectors ~f and it 

lie on the same plane. So that the relations between these frames can be given 
as follows 
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where ip is the angle between the vectors ~$ and ~N . The derivative formulae of 
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the Darboux frame is 



f 
it 





kg 

k n 



>9 





1 



where fc 9 , fc n and T ff are called the geodesic curvature, the normal curvature 
and the geodesic torsions, respectively. Here and in the following, we use "dot" 
to denote the derivative with respect to the arc length parameter of a curve. 
The relations between the geodesic curvature, normal curvature, geodesic tor- 
sion and k, t are given as follows, [11] 



kg = k cos ip , 



k n = K sin ip , Tg = T + |f 



Furthermore, the geodesic curvature k g and geodesic torsion r g of curve a (s) 
can be calculated as follows, [IT] 
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In the differential geometry of surfaces, for a curve a (s) lying on a surface S 
the following relationships are well-known, [TT] 

i- a (s) is a geodesic curve if and only if k g = 0, 

ii- a (s) is a asymptotic line if and only if k n = 0, 
hi- a (s) is a principal line if and only if r g = 0. 

Through the every point of the surface a geodesic passes in every direction. A 
geodesic is uniquely determined by an initial point and tangent at that point. 
All straight lines on a surface are geodesies. 

Along all curved geodesies the principal normal coincides with the surface nor- 
mal. Along asymptotic lines osculating planes and tangent planes coincide, 
along geodesies they are normal. Through a point of a non-developable surface 
pass two asymptotic lines which can be real or imaginary. 



3 Inextensible Flows of Curve Lying on Oriented 
Surface 

Throughout this paper, we suppose that 

a : [0,1] x [0,w) -> M C E 3 

is a one parameter family of differentiable curves on orientable surface M in E 3 , 
where I is the arclength of the initial curve. Let u be the curve parameterization 

variable, < u < I. If the speed of curve a is denoted by v 
arclength of a is 



()7t 



8 a 



then the 
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s{u) = I\\^\\ du = I 



v du. 



(3.1) 



The operator is given in terms of u by 



JL = IJL 

ds v du ' 



(3.2) 



Thus, the arclength is ds = v du. 



Definition 3.1 Let M be an orientable surface and a be a differ entiable curve 
on M in E 3 . Any flow of the curve a with respect to Darboux frame j 7 ?, if, rt| 
can be expressed following form: 



(3.3) 



Here, /i, / 2 and fz are scalar speed of the curve a. Let the arclength variation 
be 



S (w,t) = / vdu. 
o 



(3.4) 



In the Euclidean space the requirement that the curve not be subject to any 
elongation or compression can be expressed by the condition 



e-S(u,t) = J 
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at ua 



0,tie [0,1]. 



(3.5) 
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Definition 3.2 A curve evolution a(u,t) and its flow on the oriented sur- 
face M in E 3 are said to be inextensible if 



d_ 

dt 



d-ct 



du 



0. 



Now, we research the necessary and sufficient condition for inelastic curve flow. 
For this reason, we need to the following Lemma. 

Lemma 3.1 In E 3 , let M be an orientable surface and | ~T , ~g*, it | be a Dar- 
boux frame of a on M. There exists following relation between the scalar speed 
functions f\ , fi , fz and the normal curvature k n , geodesic curvature k g of a the 
curve 

% = %t~ hvk g - hvk n . (3.6) 
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Proof. Since and Mr commute and v 2 — ( C^, ) , we have 

du at \ du 7 du I 

On,&E JL / d~ct d~a 

* u dt ~ dt \ du' du 



Thus, we reach 

dv dfi 
dt du 

If we take in the conditions of being geodesic and asymptotic of a curve and 
Lemma 3.1, we give the following ■ 

Corollary 3.1 If the curve is a geodesic curve or asymptotic curve, then there 
is following eguations 

If = ^it ~ h v k n 

or 

respectively. 

Theorem 3.1 Let |t,^",^"| be Darboux frame of the curve a on M and 

= fi 7^ + f- 2 if + f^lt be a differentiahle flow of a in M 3 . Then the flow is 
inextensible if and only if 

Proof. Suppose that the curve flow is inextensible. From equations (|3.4I) and 
for u € [0, 1] , we see that 



u u 

J^S (u, t) = J ^-du = J - f 2 vk g - f 3 vk n ^j du = 0. 



(3.8) 



^ = hvk g + hvk n . (3.9) 



o 

Thus, it can be seen that 



Considering the last equation and (|3 . 2[) . we reach 

§ t=hkg + hk n . 

Conversely, following similar way as above, the proof is completed. 
From Theorem 3.1, we have following corollary. 
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Corollary 3.2 i- Let the curve a is a geodesic curve on M. Then the curve 
flow is inextensible if and only if = f3k n . 

ii- Let the curve a is a asymptotic line on M. Then the curve flow is inex- 
tensible if and only if = / 2 k g . 

Now, we restrict ourselves to the arclength parameterized curves. That is, v = 1 
and the local coordinate u corresponds to the curve arclength s. We require the 
following Lemma. 



Lemma 3.2 Let M be an orientable surface in E 3 and be a Dar- 

boux frame of the curve a on M. Then, the differentiations of | T , ~$ , ~rt\ with 
respect to t is 



where ip = ^^f-, 
Proof. Since -ii and are commutative, it seen that 

dt ds 1 

Substituting the equation ()3.7j) into the last equation and using Theorem 13.11 
we have 



Now, let us differentiate the Darboux frame with respect to t as follows; 
= -<7\f) = ( — ,?) + ( T 



dt\ * I \ dt ' 3 / \ ' 9i 
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From (|3.10p and p. lip , we have obtain 
d ~$ ft i + d 

aT =-U^ + - ea 



"^r- = - I /lfcg + — - /3T ff J T +1pn 



and 



respectively, where ip = / ^ , 
If we take into consideration last Lemma, we have following corollary. 

Corollary 3.3 Let M be an orientable surface in E 3 . 
i- If the curve a is a geodesic curve, then 

# = - (flkn + ^t + hT B ) f ~ 4>t, 

where ip = ( -S-, n 



dt 

ii- If the curve a is a asymptotic line, then 



+ f2Tg) 2* - j>t, 



dlt _ I df s 
dt \ ds 



where tp — (737-, ~f% 



Hi- If the curve is a curvature line, then 
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where ip = ^ , it 

Theorem 3.2 Suppose that the curve flow = f 1 T + f^ + f^lt is in- 
extensible on the orientable surface on M . In this case, the following partial 
differential eguation are held: 

dfcg _ dfi, , f dkg_ d 2 f 2 _ dh„ _ f drg_ _ f , _ dh _ f 2 
dt ~ ds "^"Wi ds ~ ds 2 ds '9 ■> 3 ds J lAn ' g 9s >g 32' g > 
dk n _ dhu i f dk n , d 2 h ,df 2 , fdr g , , dh f 2 

^S- — — f 1 h h - ^3. h j-f Q l. T j_M 



#n = (-flkn - ^ - /2T 3 
^ S = ("MP " H + /3T 9 ) 



y ■ 
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Proof. Since we get 



a d~f 

3s dt 



a 



( 



fit ^" + /l ^a's" + 3s2 3 + ot T 9 + / 2 
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at ds dt v9 y " " >~ dt y " 9 ot ' 

Thus, from the both of above two equations, we reach 



dt 



n — > 
n 



d n 

"~dT 



dkg_ _ dfi dkg_ 
dt - ds g + h ds ' 



d 2 f 2 df; 



dr, 



dh 



ds 2 ds 
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and 

dkn = djj_ k + , 

dt ds n ds ds 2 ds 9 ' ds 

Noting that = Jj^f-, it is seen that 
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(3.15) 



From the equations (|3.12p and (|3.14j) . it is seen that 

= - fik n T & 

= (-flkn 



#n = -fik n T g - ^fr g - f 2 r g 2 
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By same way as above and considering J| = , we reach 



dk n dfx dk n d 2 h . 9f 2 , , dr g , ( 

Hence, from the equations (|3.13p and (|3. 161) . we get 



Qf 2 



Thus, we give the following corollary from last theorem. 

Corollary 3.4 Let M be an orientable surface in E 3 . 

i- If the curve a is a geodesic curve on M, then we have 



dk n dfx dk n d 2 f 3 df 2 . dr g df 2 , 2 

~dt = ~~ds ~ds~ ~ds 2 ~ ~ds~ Ta ~ds ~ds~ Ta ~ 3 



dr g df 2 dtp 

-77- = TT- K + hk n T g + — 

at os as 



ipk n = (-fik n - ^ - f 2 T 6 
ii- If the curve a is a asymptotic line, we have 



dk g dfi dk g d 2 f 2 df 3 dr g df 3 2 

~~dt ~ ~d7 9 ^~d7 ~ol 2 ~ ~ ~ds~ Tg ~ ^~ds~ ~ ~ds~ Ta ~ ^ 

dr g dtp 
~dt^~ds~ 

and 

Ipkg = (-flkg - ^ + fiTg ^ Tg. 

Hi- If the curve a is a curvature line, then we have 

9kg _ 8/1 u , f dkg_ d^h 

dt ~ ds 9 -lids "T" df 

dk n _ dfi r. , f dk n i 8 h 

dt ~ ds ' J 1 Qs T ds 2 ■ 
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